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Abstract 

We propose to analyse the statistical properties of a sequence of 
vectors using the spectrum of the associated Gram matrix. Such se- 
quences arise e.g. by the repeated action of a deterministic kicked 
quantum dynamics on an initial condition or by a random process. 
We argue that, when the number of time-steps, suitably scaled with 
respect to h, increases, the limiting eigenvalue distribution of the Gram 
matrix reflects the possible quantum chaoticity of the original system 
as it tends to its classical limit. This idea is subsequently applied to 
study the long-time properties of sequences of random vectors at the 
time scale of the dimension of the Hilbert space of available states. 
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Discretising a classical dynamical system is in order if we want to simulate 
it on a computer. Its compact phase space may for this purpose be covered 
by a large number of small patches of Lebesgue measure The evo- 

lution, which we assume measure preserving and discrete in time, translates 
approximately in a bijection of the patches. Such a description involves al- 
ways an approximation as patches change their shape in the course of time. 
It is of course helpful for actual model systems to label the patches in a way 
that mimics the kinematic structure of phase space. So, we obtain after this 
coarse-graining procedure for each a one to one transformation tt of the set 
{1,2,... , N} that determines the evolution during one tick of the clock. The 
phase portrait consists in partitioning the discrete phase space into closed 
orbits of vr and the crucial information is the number of orbits together with 
their length as a function of N. An ergodic island of non-zero measure in 
the dynamical system will signal its presence by the occurrence of an orbit 
with a period proportional to A^. Iterating the dynamical map on an initial 

point io provides us with a sequence z = ^io, 7r(io), 7r^(io), • • • ^ of points in 
{1,2,... , N} and we can distinguish between points belonging to ergodic or 
regular regions of phase space by examining the period of the time sequence 
of io as a function of N. 

Truly quantum dynamical systems with compact phase space are finite di- 
mensional in virtue of the uncertainty principle. As each state occupies a 
same volume h, the dimension of their Hilbert space of states is 1/h. Planck's 
constant has here a rather symbolic meaning: for rf-dimensional systems it 
is the d-th power of the actual Planck constant, while h = l/(2j -|- 1) for a 
spin with angular momentum j. The d-dimensional complex Hilbert space, 
or, more precisely, the space of complex rays in is the quantum space 
with d elements. The space of rays is called the projective Hilbert space of 
dimension d and denoted by prC^. In dimension 2, it turns out to be the 
unit sphere in with antipodal points identified. In contrast to a classical 
space, distinct points can lie arbitrarily close, the distance between the rays 
[</?] = C(fi and — Cip generated by the normalised vectors (p and ip being 

d(M,[V']):= inf \\<p - z ^\\ ^ 2 - 2\{^,^^;)\ ^ Asm' ^ , (1) 

\ / z€C, \z\=l I 

with Q G [0,7r/2] the angle between the rays. The maximal separation be- 
tween points is reached when they correspond to orthogonal rays. Projective 
Hilbert spaces carry a natural Riemannian structure given by the Study- 
Fubini metric, but we shall not be so much concerned here with this contin- 
uum feature and rather focus on their discrete aspects. 
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A quantum evolution in discrete time, also called kicked evolution, is deter- 
mined by a unitary Floquet operator u. In Schrodinger picture, ^-^ uip is 
the evolution between two consecutive kicks. We use the same notation to 
denote the corresponding evolution in the space of rays: p \—>- up. We now 

face the problem of studying time sequences p = (^po, upo, u'^ po, . . gen- 
erated by a Floquet operator u as it acts repeatedly on an initial condition 
Po- 

In the vast literature on quantum chaos, dynamical properties of quantum 
systems are often investigated by considering the temporal behaviour of the 
Husimi or the Wigner functions corresponding to well-localised states in 
phase space, such as coherent states. Such a description relies, in taking the 
(semi) -classical limit, on a definite geometrical picture of the corresponding 
classical phase space i.e. on a particular choice of basic observables such as 
the usual position and momentum or angular momentum. Many references 
can be found in |l| and We argue in this letter that it is worthwhile to 
put things in a more abstract perspective. Gram's matrix, or rather its spec- 
trum, provides us with a powerful tool to analyse the statistical properties 
of time-sequences of points in a projective Hilbert space. 

The Gram matrix G{ip) of a sequence (p = ^v^(l), V^(2), . . . , ^p{K)j of vectors 
is 
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G(c^) is positive definite and its rank equals the dimension of the space 
spanned by the ^{j) 0. In particular is linearly independent if and only 
if det{G{(fi)) 7^ 0. The spectrum of G(yj) is independent on the order of the 
in (p and on multiplying the ^{j) with a complex number of modulus 

1. This means that for a given sequence p = ^p(l), p(2), . . . , p{K)j of 

points in a projective Hilbert space, specified by normahsed vectors (p{j) 
as p{j) = [<p{j)], the spectrum of G{(p) depends only on p and that it is 
insensitive to the order of the points in p. It may therefore be denoted by 

Let us for a moment consider a classical word i = . . . , i{K)j 

where the letters are chosen from a given alphabet {1,2,...}. In fact, i is in 



3 



one to one correspondence with a sequence p = \ [ei(i)], [ei(2)], • • • , j 
of vectors through the identification of j with ej for an orthonormal basis 
{ei, 62, . . . } of a Hilbert space. Grouping with equal index j, the Gram 
matrix is block-diagonal with block E{j) of the type 



The dimension of E{j) is precisely the multiplicity m{j) of j in As the 
spectrum of E{j) consists of the non-degenerate eigenvalue m{j) and the 
m(j) — 1 degenerated eigenvalue 0, we find that determines precisely 

the amount of different numbers appearing in i with their multiplicity i.e. 
the relative frequencies of the different letters in i. The spectrum of the 
Gram matrix of a very regular sequence i will consist of a few large naturals 
and a highly degenerated while for sequences with many different indices 
the spectrum will be concentrated on small natural numbers appearing with 
high multiplicities. A same interpretation remains valid for the general non- 
commutative Gram matrix with spectrum concentrated around small 
natural numbers points at a vector wandering wildly through the Hilbert 
space of the system and is therefore a sign of chaotic behaviour. More regular 
motion, such as precession or slow diffusion, signals its presence by large 
eigenvalues and a high occurrence of eigenvalues close to 0. In contrast to the 
classical case however, eigenvalues are no longer limited to natural numbers 
so that a same point in a projective Hilbert space can now be visited a 
fractional number of times. 

The existence for quantum systems of an intermediate time scale with in- 
teresting and describable behaviour in between the very short one, of order 
— log h, where the quantum system slavishly follows its classical limit and the 
very long one where quasi-periodic behaviour, due to the discreteness of the 
spectrum of the Floquet operator, is dominant, is a central theme in many 
papers [Q. We are, more precisely, interested in the limiting eigenvalue dis- 
tribution of the Gram matrix when its dimension, the number of time-steps, 
tends to infinity, appropriately scaled with respect to the quantum parame- 
ter. To obtain this, we consider the limit of the empirical measure 



where the jj are the eigenvalues of the Gram matrix. The limiting distribu- 
tion should refiect the quantum character of the dynamics as it tends to its 
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classical limit. We do not claim that we can settle this question but we shall 
at least provide some rigorous support for this expectation. 

Instead of considering a genuine unitary dynamics acting on an initial con- 
dition, we consider sequences p = ^p(l), p(2), . . . , p{K)j of points in the 
projective Hilbert space of dimension A^, independently and randomly chosen 
with respect to the uniform measure. Recall that prC^ is compact and that 
it carries a unique normalised measure, called uniform, which is invariant 
under the action of every N x N unitary. Picking independent normalised 
vectors randomly with respect to this measure is quite different from picking 
the components of the vectors with respect to a given basis in an independent 
and random way with respect to some suitably chosen probability measure. 
Next, we compute the spectrum of the Gram matrix of such a sequence. This 
spectrum is of course a random object but it turns out that, in the limit of 
large N and for a rescaled time r = K/N, the spectral distribution tends to a 
definite limit given by the Marchenko-Pastur distribution fij. [Q. The actual 
computations are somewhat involved and will be presented in Though 
the Gram matrices are given in terms of independent random vectors, there 
is no independence between the entries. E.g., each Gram matrix is positive 
definite which is totally incompatible with independence of matrix elements. 
This makes the computation quite different from Wigner's random matrix 
computation, fi^ is obtained either by a combinatorial argument in terms of 
its moments or by determining the expectation of the resolvent of the Gram 
matrix. 

When < r < 1 the probability measure fir is given by a continuous density 
p. For very small r, we must choose relatively few vectors in a large space. 
This will often lead to almost orthogonal choices and therefore p will be 
concentrated around 1. When r increases to 1, there is a fair chance of many 
vectors overlapping and the support of p will simultaneously extend towards 
0, which is a lower bound of its support, and to larger positive values. When 
r > 1 there will almost surely be a sizeable degree of linear dependence 
responsible for a high multiplicity of the eigenvalue 0. In fact, it turns out 
that pr decomposes for r > 1 into an atom at and an absolutely continuous 
part: 



where p is the absolutely continuous part of Pr- The weight of p is 1/r. 
Moreover p is compactly supported in the interval [(v^" 1)^? (v^+ 1)^]- A 
similar computation in the classical case yields a Poisson distribution. This 
result is reminiscent of Wigner's semicircular distribution for the spectrum 
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of large random matrices where compactness of the hmiting distribution is 
also a typical feature of the non-commutativity . A simple measure of the 
dynamical entropy of the system is the length of the support of fi^-. For our 
random dynamics, this quantity grows as 4:^/t, in contrast to an expanding 
chaotic dynamics where the entropy grows linearly in r. 
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Figure 1: The limiting spectral distribution of G{p) 



The figure shows the limiting spectral distribution of Gram matrices in the 
region 0.02 < r < 3. The 6 contribution of weight (r — l)/r that appears for 
r > 1 is rendered by the fat line, which has height (r — l)/r. The continuous 
part is, for all values of r, only non-vanishing for {y/r — 1)^ < x < (a/t + 1)^, 
but this is only visible in the figure for moderately small values of r. For 
r tending to 0, a 5 distribution at x = 1 will appear and the probability 
density has for r = 1 a singularity at x = of order —1/2. 
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Figure Captions 

Figure 1: The hmiting spectral distribution of G{p) 
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